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Abstract 

> 

m . 

<^ , The dynamics of relativistic stars and black holes are often studied in terms 

of the quasinormal modes (QNM's) of the Klein-Gordon (KG) equation with 



different effective potentials V(x). In this paper we present a systematic study 
of the relation between the structure of the QNM's of the KG equation and the 
form of V(x). In particular, we determine the requirements on V(x) in order 



for the QNM's to form complete sets, and discuss in what sense they form 

H : 

complete sets. Among other implications, this study opens up the possibility 

of using QNM expansions to analyse the behavior of waves in relativistic 

systems, even for systems whose QNM's do not form a complete set. For such 

systems, we show that a complete set of QNM's can often be obtained by 

introducing an infinitesimal change in the effective potential. 
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Typeset using REVTgX 



1 



I. INTRODUCTION 



Gravitational waves have been of theoretical interest since the appearance of the general 
theory of relativity, and of experimental interest for several decades. The generation and 
propagation of gravitational waves are often discussed in terms of a source region, a strong- 
field and a weak- field near zone, and a local and distant wave zone [[1]]. For the purpose 
of this paper, a simpler point of view will be useful for understanding: (i) a source region 
where the radiation is generated, (ii) an intervening region with a non-trivial geometry of 
the background spacetime, in which the propagation of waves is significantly different from 
that in flat space, and (iii) an "outside" region where the spacetime can be taken to be 
flat, in which the observer is located. A close analogy is a molecule or molecules [region (i)] 
radiating in an optical cavity [region (ii)]. The cavity significantly modifies the character of 
the electromagnetic radiation observed outside [region (iii)]. 

The interest of this paper is on the effects of the intervening nontrivial spacetime, i.e., 
the "cavity", in gravitational systems. For example, consider radiation from a compact 
astrophysical configuration, e.g., a binary neutron star system, embedded in a galaxy. The 
source region is the binary system, the "cavity" is the galaxy, and the rest the "outside". 
The presence of a "cavity" in principle could have a significant influence on the radiation 
observed from the outside. The issue of a recent debate [0] is whether the "cavity" , in this 
case, the galaxy, can suppress, by many orders of magnitude, the emergence of radiation 
generated by the source ||. Another example is a particle radially plunging into a black 
hole of mass M. In this case the "source region" and the "cavity" have the same size (as is 
true for a laser cavity). The "source" is significant only within ~ 10M of the plunge, and 
the effective potential (in the perturbation theory of a Schwarzschild hole [|J ) representing 
the "cavity" is also significant only within ~ 10M. It is well-known that in this case the 
"cavity" modifies the observed radiation significantly. 

Radiation from optical cavities is often discussed in terms of the "modes" of the cavity. 
Because the waves escape to infinity, the total energy in the finite part of the system decreases 
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with time. Thus, these "modes" are characterized by complex frequencies ujj, with Im uo,j < 0. 
The corresponding eigenfunctions are defined by the outgoing wave condition at infinity. 
These "modes" are in fact quasinormal modes (QNM's). The observation of the QNM's 
from the outside gives directly information on the spatial structure of the cavity (but not 
the structure of the source, unless a particular model is assumed); e.g., in obvious notation, 
Uj ~ jnc/L, where L is the length of a simple 1-dimensional optical cavity. 

The same discussion extends to the gravitational case. The QNM's of black holes @] 
and relativistic stars f| have been subjects of much study. In numerical simulations, it is 
often found that the radiation observed in many black hole processes is dominated by the 
QNM's of the background spacetime associated with the hole. In view of the optical analogy, 
this is not surprising — the distant observer sees only the QNM's of the laser cavity, but 
not the details of the source. The fact that gravitational QNM's tend to be more "leaky" 
does not change this qualitative understanding. 

It is therefore particularly exciting that gravitational QNM radiation may be observed by 
LIGO, VIRGO, and other detectors |7| in the next decade. These observations will provide 
the possibility of "seeing" directly interesting spacetime structures of various gravitational 
systems, e.g., the spacetime around a black hole — in much the same way as the spectrum 
of a laser permits a distant observer to infer the geometry of the cavity. 

Although the QNM's of black holes are quite well studied, the general behavior of QNM's 
in more complicated gravitational systems has not received much attention. For example, 
in view of the possible detection of gravitational waves, one would like to know: How much 
would the QNM frequencies of a black hole be shifted by the astrophysical environment that 
it is in, e.g., an accretion disk or the host galaxy? Can one calculate the shift perturbatively? 
How much would each QNM be excited in an astrophysical process, e.g., a star plunging 
into a massive hole? There are many such questions of both theoretical and observational 
interest . 

In an ongoing project, we study the general properties of QNM's in gravitational systems. 
In a recent letter |§, we outlined the results obtained in studying the QNM's of the Klein- 
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Gordon (KG) equation 




dx 2 
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+ V(x) <f)(x,t) = 0, 



(1.1) 



which is often used to describe wave propagation in curved space, with V(x) describing the 
scattering of waves by the geometry. In this paper, we provide the details of the results 
reported in |§. We also present other related results and discuss the implications. 

For a set of modes, the first question is whether they are complete. Completeness can 
mean several different things, and we shall discuss the relations among the various senses of 
completeness in Sec. Ill; for the moment we take the simple point of view that the QNM's 
are complete if the evolution of the wavefunction can be expressed, in some domain of time 
and space, as a sum 
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where fj(x) are eigenfunctions with complex frequencies Uj. In particular, if this holds at 
t = 0, then an arbitrary function of x within a certain domain should be expressible as 



Both of these properties can be discussed in terms of a Green's function G(x, y; t) for propa- 
gating a source at y to an observation point at x in time t. Our purpose is to elucidate under 
what conditions would completeness hold. This problem is nontrivial because usual proofs 
of completeness for the normal modes of conservative systems rely crucially on hermiticity. 
It is therefore particularly interesting that Price and Husain || were able to give a model of 
relativistic stellar oscillations in which the QNM's do form a complete set. It was suggested 
H that the completeness of the QNM's stems either from the equality of the damping times 
of all the QNM's, or from the absence of dispersion and backscatter in that model. We shall 
see in this paper that these properties are not the controlling factors of completeness. 

There is a limit in which QNM's should behave trivially, like the normal modes of her- 
mitian systems, namely the limit in which the leakage goes to zero; mathematically this 
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means that the quality factor Qj ~ Re ouj/(2 |Im ujj |) — > oo. In this limit, the QNM's 
become normal modes and are complete. In the present paper we study the completeness 
of gravitational QNM's whose quality factors are not large, and we are interested in results 
that survive to all orders in the leakage, which may be characterized by 1/Qj- We show that 
the QNM's are complete in a broad class of models quite independent of leakage, provided 
the following three conditions hold: 

(i) The effective potential V(x) in (1.1) is everywhere finite, and vanishes sufficiently 
rapidly, in a sense to be defined in Sec. II, as \x\ — > oo; 

(ii) there are spatial discontinuities in V(x) demarcating the boundaries of a "cavity"; and 

(iii) consideration is limited to certain domains of space and time, which will be spelt out 
below. 

It is useful to comment heuristically on these conditions at the outset. If (i) does not hold, 
i.e., if V(x) has a significant tail at large x, then it is possible for a disturbance originating 
at y to propagate to a large x' 3> x, and be scattered back to the observation point x. The 
time taken would be 0(x'), and since x' can be arbitrarily large, this could lead to a late time 
tail in the wavefunction; crudely speaking this would be t~ a if asymptotically V(x) ~ x~ a 



(though this is modified in the presence of a centrifugal barrier [Kj). Such a late time tail 
would not be described by a sum of discrete QNM's, which must behave exponentially at 
large times. Thus, the QNM's cannot be complete for V(x) having a tail. The late time 
tail, and its relationship with the spatial asymptotics of V(x), is now thoroughly understood 

A discrete sum over QNM's (like a Fourier series rather than a Fourier integral) cannot 
possibly be complete over all space; completeness, if it holds at all, can only be limited to 
the inside of a finite interval. The need for a natural demarcation of this interval "explains" 
why discontinuities are needed. In the model of M there is such a discontinuity in the speed 
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of wave propagation at the model stellar surface. This also explains why the completeness 
relation (1.3) must be limited spatially. 

The need to also limit the time domain in some cases can be understood from the 
following example. Suppose the observation point is far from the source. Then there will 
be transients propagating directly from the source point y to the observation point x (i.e., 
the high frequency components), without being much affected by the potential; these have 
nothing to do with the QNM's. Only past a certain finite prompt time t p = 0(x), when 
these transients have passed, would the QNM's be complete. This case is most relevant for 
the discussion of gravitational waves seen by a distant observer, and numerical experiments 
indeed show that quasinormal ringing dominates only after the passage of transients. 

It should be emphasized that the broad class of models satisfying the above conditions 
and hence exhibiting completeness allow dispersion, backscatter, as well as differences in the 
damping times of the QNM's, showing that these are not essential ingredients. This broad 
class of models also allow the leakage to be not small; in other words, our results are valid 
to all orders in l/Qj. 

There is one final technical hedge. The QNM sums are in fact often divergent series, 
but converge to the correct answer if regularized in a standard way. The simplest statement 
of regulation is: keep only modes with Re Uj > by symmetry argument, then replace all 
times t by t — ir, with r — > + . We shall come to the need for this in Sec. II, and we 
shall also demonstrate, through a numerical example in Sec. IV, that the regulated series is 
useful in practice. 

The general need for regularization is a major difference between the completeness prob- 
lem of the KG equation (1.1) and the wave equation 

d 2 d 2 



n 2 (z) 



dt 2 dz 2 



r(>(z,t) = 0. (1.4) 



where n(z) may be regarded as a position-dependent refractive index in a scalar model 
of electromagnetism. We have previously studied the completeness of QNM's in the wave 
equation case [|TT|1 . In f9j, the stellar oscillation model is described by (1.4), with n(z) = 1 for 
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< z < L representing the interior of the star, and n(z) a non-zero constant for L < z < oo 
representing space exterior to the star. The boundary conditions are ip(z = 0) = and 
outgoing waves at infinity. 

In this paper we work in terms of the KG equation (1.1) since it is more realistic as a 
model for wave propagation in curved space, and has been much studied for this purpose. 
For example, in a static, spherically symmetric spacetime, 

ds 2 = g tt {r)dt 2 + g„(r)dr 2 + r 2 (d6 2 + sin 2 6dip 2 ), (1.5) 

a KG scalar field $ can be expressed as 

9 = Y t -<f>im{x,t)Y lm {d,ip)- (1-6) 

Im r 

The evolution of each (/>i m (x,t) is given by (1.1) with the effective potential 

yf x \ = g JV + _ J_J?E(J_^ _ _]_^9rr^^ /j ^ 

r 2 2rg rr g tt dr g rr dr 
The variable x in the KG equation is related to the circumferential radius r by 



x = [J-tLdr. (1.8) 
J V 9tt 

The particular case of (1.1) arising from a Schwarzschild spacetime is well studied. The 
Maxwell field and the linearized gravitational waves satisfy the same equation with slightly 
different potentials V(x) [§]. 

The wave equation (1.4) is related to the KG equation (1.1) by a transformation: 

dx/dz= n(z), 

iP= n" 1/2 0, (1.9) 

with the potential related to n(z) by 

V = (2n 3 )- 1 (d 2 n/dz 2 ) - (3/4n 4 )(dn/dz) 2 . (1.10) 

Hence the results obtained in one case can be restated in the other, up to possible compli- 
cation coming from regularization, which is investigated in this paper. As we shall see later, 



discontinuity is a crucial feature in the consideration. We note that V(x) with a discontinuity 
in its p-th derivative maps into n(z) with a discontinuity in its (p + 2)-th derivative. 

Our results for discontinuous potentials, in addition to clarifying the conditions for com- 
pleteness, are significant in two ways. First, they are directly applicable to physical models 
with a discontinuity, e.g., a stellar surface as in the model of ||. Second, any smooth po- 
tential can be approximated, to arbitrary precision, by a discontinuous one. In fact, any 
numerical scheme that employs finite spatial differencing has in effect replaced V(x) by one 
with a discontinuity in a high order derivative. Thus the complete set of QNM's of the 
discontinuous system can be used as an effective calculational tool for the original system. 
Further remarks on this aspect are given in Sec. V, while detailed studies along this line 
will be reported elsewhere. 

In Sec. II we develop the formalism via the Green's function in the complex frequency 
plane, and demonstrate that for a class of discontinuous potentials V, the Green's function 
is given by a sum of QNM's. The conceptual analogy with optical cavities extends to a close 
parallel in the mathematical treatment with the optical case but there are technical 
differences which are stated in this paper. Sec. Ill spells out more clearly the different 
senses of completeness and also discusses briefly the extension to KG equation with a non- 
zero mass. In Sec. IV, we verify numerically the completeness relation in a model problem. 
The paper ends with some discussions and conclusion. 
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II. THE GREEN'S FUNCTION 



A. Green's function in frequency plane 

We model the propagation of waves in curved space by the KG equation (1.1). The 
potential V(x) is assumed to be bounded and positive, and represents the scattering of waves 
by the background. Generically, we take V a = lim^oo V(x) = 0. The case V a 7^ 0, which 
represents a massive field, will be discussed briefly in Sec. III. The spatial coordinate x often 
represents a radial variable [cf. (1.8)], so we shall first consider a half-line problem (x > 0) 
with the regular boundary condition <f>{x = 0, t) = 0. The full line problem (— oo < x < oo), 
e.g., x being the tortoise coordinate in the Schwarzschild case, will be studied later. The 
QNM's are eigen-solutions 4>(x,t) = <p(x)e~' luJt of (1.1), where satisfies 

D$ = i-u 2 + V(xMx) = 0, (2.1) 

with <p(x = 0) = and the outgoing wave boundary condition at infinity. The eigenf unctions 
and eigenvalues are defined as 4>{x) = fj(x) and to = Uj. 

The Green's function G(x, y; t) for the system is defined by DG = 8{x — y)5(t) with the 
initial condition G = for t < and the boundary conditions that (i) G — for either 
x = or y = 0, and (ii) the outgoing wave condition as either x — > oo or y — > oo. The 
corresponding Green's function in the frequency domain is 

DG(x,y;uj) = 6(x-y). (2.2) 

The strategy is to express G in terms of G, and attempt to close the contour by a large 
semicircle in the lower half u plane. 

Introduce two auxiliary functions f(u,x) and g(u,x), which are the solutions to the 
homogeneous time-independent KG equation Df(u,x) = Dg(u>,x) = with the boundary 
conditions f(u,x = 0) = 0; f'(uj,x = 0) = 1 [12|] and lim^oo [g(uJ, x) exp(— iujx)] = 1. With 
these auxiliary functions the Green's function is then given by 
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G(x,y;u) = f(w,x)g(u,y)/W(u) for < x < y, 

= f{u,y)g{u,x)/W{u) for < y < x, (2.3) 

where the Wronskian W(u) = g(u,x)f'(u,x) — f(uj,x)g'(uj,x) (' = d/dx) is independent of 
x. G(x,y;ui), of course, may be singular in u at the singularities of / and g; otherwise, it 
is analytic except at zeros of W(u). At these zeros / and g are proportional to each other, 
so each of them satisfies the regular boundary condition at x — 0, and the outgoing wave 
boundary condition as x — > oo. Hence these frequencies are exactly the QNM frequencies 
ujj and f(ujj,x) = fj(x). We shall further assume for simplicity that these zeros are simple 
|13| , so that the residues of G(x, y; uj) at these zeros are given by 

Kj = ffa, x)g{u j , y)/[dW{u = u^/du). (2.4) 

Multiple zeros can be handled readily. To study the physical meaning of the denominator 
dW/duj, start with the defining equation for f(uj,x) and g(u,x). The usual manipulations 
lead to 

f x 

(u 2 -u]) dx f(uj,x)g(u,x) = g(u,x)f'(uj,x) - gf(u,x)f(uj,x)\ Q , (2.5) 

where the integral is taken along any contour from x = to x = X. Since both f(uj, x) and 
g(u,x) are outgoing waves at x — X, the right-hand side of (2.5) becomes 

i(oOj - u)f{ujj, X)g(u, X) - [g(u 7 0)f(uj, 0) - g'(oj, 0)f(uj, 0)] , (2.6) 

as X — > oo. Differentiating (2.5) with respect to u and taking uj — > ujj then gives 

&J f*«fc f^zH^z) + ±-K^X) 9( ^X) = -±. dW ^ g = ^ . (2.7) 

Since / and g are proportional to each other at these poles, the residues can be expressed 
in terms of a generalized norm of the QNM's 

< fj I fi » = l im I"" dx /(^i, x)f{uj, x) + ^-fiujj, X)f{uj, X) (2.8) 

as follows 
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K IMIM (29) 

K >~ 2u, «f J \f J » {2J) 



The norm (2.8) has been introduced in other contexts [Tl|JT^|T^] and has the following 
significant properties: 

(i) it invovles f 2 rather than \f\ 2 , and is therefore in general complex; 

(ii) there is a surface term; 

(iii) although each of the two terms on the right of (2.8) does not have a limit, the limit 
exists for the sum; and 

(iv) it reduces to the usual norm if the leakage goes to zero. 

The property (iii) is easily verified by differentiating the right hand side of (2.8) with 
respect to X, and using the outgoing wave condition, which also applies to / at u — Uj. 
The appearance of the surface term is closely connected with the fact that the momentum 
operator is not hermitian. While (2.8), with the integral taken along the real axis, is often 
convenient for actual evaluation, we may choose a deformed contour L running from x = 
to x = — oo. The first part L\ = (0, a), and the second part L<i is shown in Fig. 1. On L2, 
f(ujj,x) is defined as the outside solution (which is analytic in x for x > a), analytically 
continued to the complex plane. Along this contour, f(uij,x) — > as x — ■> —00, and we get, 
more compactly 

< fj I > = dxf(uj, x fdx } (2.10) 

making the formal analogy to the hermitian case more transparent. However, analytic 
continuation out into the complex x plane is not strictly necessary if the surface term is 
retained, as in (2.8). 

Next write G in terms of G; then upon closure of the contour in the lower half u plane, 
one sees that 

- j U 3 <<; - J 3 I J 3 > 
11 



In (2.11), the sum comes from the zeros of the Wronskian W(u) inside the semicircle, J c is 
the integral along a semicircle at infinity, and I s comes from the singularities of f(u,x) and 
g(uj,x) (see Fig. 2). The crux of the proof of completeness then lies in (i) determining the 
conditions on V(x) under which / and g have no singularities in the to plane, in which case 
I s = 0; and (ii) showing (in the next subsection) that I c vanishes if there is a discontinuity 
in the potential V(x) at some x = a > 0. Under these conditions, the Green's function can 
be represented exactly in terms of the QNM's, which then establishes these QNM's as a 
sufficient basis for discussing the dynamics. The behavior is a discrete sum of exponentials, 
and in particular, there would be no power-law dependence at large times. 

Since QNM's appear in pairs at oj — Uj and u = — u;J, it is convenient to rewrite the 
QNM sum (when I s = I c = 0) as 



G(x, y; t) = Re 



/~o - « I j I j > 



(2.12) 



where the notation j > is a shorthand for Re Uj > 0. This latter form is somewhat 
more convenient for regularization when the sum diverges. The corresponding statements 
for G(x, y; t), and consequently a QNM representation of 8(x — y) = G(x, y;t = + ), follow 
relatively simply and will be discussed in Sec. III. 

Now to determine the conditions under which f(u,x) and g(u,x) have no singularities 
in w, we appeal to well known results in the quantum theory of scattering. The defining 
equation for / and g is identical to the Schrodinger equation for a particle with mass = 1/2 
and energy = uo 2 moving in a potential well V(x). It has been proved jl6) that / and g are 



analytic functions of u if the potential is bounded and "has no tail" , in the sense that 

/ dx x I V(x) |< oo, (2.13a) 
Jo 

roo 

/ dx x e ax | V(x) |< oo for any a > 0. (2.13b) 
Jo 

Note that if condition (2.13b) is violated for some a > a a > 0, then g(u,x) may not be 
analytic for Im uo < —a Q . The singularities that appear for potentials that have a tail, 
e.g., inverse-power-law potentials, have been studied, with a focus on understanding the 
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associated late time behavior Jl0| . 

It then remains to investigate the behavior of G on the large semicircle. 



B. Asymptotic behavior of Green's function 

We first give a simple (and not strictly rigorous) derivation to highlight the essential 
ideas. The asymptotic behavior of G at high frequencies can be determined using the WKB 
method. Let <fi{x) = exp[iS , (x)] be a solution of the time-independent KG equation; then in 
lowest order approximation 

S{x) = ± r k(x')dx', (2.14) 



where the position-dependent wave number k(x) is given by k(x) = [u 2 - V(x)} 1 / 2 . The two 
auxiliary functions / and g, and hence the Green's function G, can be obtained in terms of 
4>{x) and k(x). 

Consider a potential with a discontinuity at x — a. In this situation the WKB ap- 
proximation breaks down at the discontinuity; however, one can join the two approximate 
solutions across the discontinuity using the reflection coefficient 

g'(a-) - 

H ~ [S'(a-)]* + S'{a+) ■ { } 

The discontinuity can be in any derivative of V, and R behaves as some inverse power of 
uj as \uj\ — > oo. (This will be spelt out precisely below.) We divide the discussion into two 
cases: (a) < y < a < x\ (b) < y < x < a. Referring to the discontinuity at x = a as the 
stellar surface for simplicity, we may say that case (a) is relevant when gravitational waves 
from a source inside the star is detected by a distant observer outside the star, while case 
(b) is relevant for discussing the internal dynamics of the star. 

For ), it is straightforward to show that 

G(x, y; u) ~ . v ; — 1 v n , y < a < x, (2.16) 

^k{x)k{y)[e- u ^ a ) + Re u ^ a )} 
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where 

I(u,v) — I k(x)dx « u(v — u). (2-17) 



Now consider G on the semicircle uj = Ur + iioj = Ce ,tt < < 2ir, with C — ► oo. As 
cjj — > — oo, both the numerator and the denominator of G are dominated by the term 
proportional to R, and 

GfoyjuOe-*" ~ 1 + J ^ • (2.18) 

As C — > oo, this vanishes for t > t p = max(a; + y — 2a, 0), since i? varies as an inverse power 
of u. Thus we have proved (modulo a technical complication to be mended below), that for 
a discontinuous potential, completeness in the sense (1.2) holds for t > t p (x,y). 
Next consider case (b). Instead of (2.16), one now has 

G(x y co) * MmvW-* 1 ™ + R* aM ] (2 ig) 

^/kixjk^le-^a) +Re ima)] 

Again, both the numerator and the denominator in G are dominated by the term propor- 
tional to R, and 

e -iu(t+x-y) 

G(x, y; u)e~ iuJt ~ , y < x < a. (2.20) 

uj 

As C —>■ oo, this vanishes for all t > 0. Thus completeness is again proved. 

The above derivation shows the essence of the proof, which relies on the vanishing of an 
exponential factor e _ ' w/ '°" where a > (say a = t — x — y + 2a), when to = Ce ie , n < 9 < 2ir 
and C — > oo. This derivation however must be mended in the domain | uj \— O(logC), i.e., 
within an angle A# ~ log C/C of the real axis, in which t~\ Wl \ a behaves only as a power C°~, 
rather than as an exponential []T7|. The problem is closely related to the QNM's, which lie 



precisely in this domain. 

To discuss this in a somewhat more general context, consider a potential with a disconti- 
nuity in d p V(x)/dx p at x = a, for some p > 0. Then as u — > oo, R ~ Au~ q , where q = p + 2. 
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We first determine the asymptotic position of the QNM's, which are given by the zeros of 
the denominator G in (2.16): 

e -il(0,a) + Re il(0,a) _ e -iua + J^-q^a = 0> ( 2 .21) 

with the solution 

u ja — jit — i [q log(jvr) - log{-a q A)} /2, (2.22) 

where j takes an integer values. The distribution of QNM's is shown schematically in Fig. 
3. 

The derivation sketched above relies on the fact that the integrand vanishes on the large 
semicircle as C — > oo. However, for u ~ Uj, one can readily obtain the following estimate 
for the integrand: 

\T3\{t-x-y)/2a I -q(t-x-y)/2a\ 

\G(x, y; u)e-^\ ~ l -^ T ~ l - n l - , (2.23) 

\oj\ \u>\ 

for both ) and case (b). In order to justify the derivation for the completeness 

of QNM, one could restrict t > x + y, but this would greatly limit the usefulness of the 
completeness relation. However, it will be shown in the next subsection that this difficulty 
can be surmounted by using a standard regularization scheme. 

The situation for t < x + y in case (b), which does not give convergence, is markedly 
different from the case of the wave equation with a step discontinuity, which was considered 



by us in the context of optics [0], and also by Price and Hunain as a model for gravitational 
waves 0. In that case, the analog of the potential V(x) is u 2 e(x), where e(x) = n(x) 2 
is the dielectric constant. Any discontinuity in e(x) is therefore amplified as uj — > oo by 
two powers; specifically, if there is a discontinuity in d p e(x)/dx p , then R ~ Au~ q at high 
frequencies, with q = p. The difference by two powers also follows from the mapping (1.10). 
Thus for p = 0, G is sufficiently bounded for the proof of completeness to go through fllTf , 
and the statements about the possibility of expressing G and dG/dt as sums over QNM's 
are strictly valid for the wave equation [jlBf . If p > 1, then the same problem arises even for 
the wave equation. 
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However, both the case of the KG equation, and the case of the wave equation with 
p > 1, can be rectified by a simple scheme of regularization. 

C. Regularization 

First, it follows directly from the reflection symmetry relation G(x, y\ — oj*) = G(x, y\ u)* 
that 

roo Aijj 

G(x,y;t) = 2Re —e^Gix^u). (2.24) 

JO 2,71 

As the integrand is well behaved at both endpoints, we can multiply it by an extra factor 
e _wr , r > 0, and then take the limit r — > + afterwards, i.e., 

G(x, y; t) = 2 Re lim / — e~ iw{t ~ iT) G(x, y; to). (2.25) 
r^o+ Jo 2ix 

This is equivalent to assigning t an infinitesimal negative imaginary part, and provides 
sufficient regulation on the part of the contour with \ujj\ small to make the integral vanish 
as the contour expands to infinity in the lower half plane. 

Second, we deform the integration contour in (2.25) into a sequence of expanding rectan- 
gular contours Y n = T nl \J T n2 U r n3 in the lower half uj plane as shown in Fig. 4. However, 
instead of letting the width of the contour T n expand continuously to infinity (which would 
hit the QNM's), each Y n is defined to cut roughly mid- way between ui n and ui n+ i, and even- 
tually we take the limit n —>■ oo. The height Q n of the contour is chosen to be proportional 
to n. As a consequence, G can be expressed as a sum over QNM's lying within the rectangle 
formed by the real axis and T n |T{|, plus an integral along the contour T n : 

fjU)f ,(!/)< '~ J : ' 



G(x, y; t) = 2 Re / — e"^ ( *- ir) G(x, y; to) + Re lim 
jr„ 2tc t^o+ 



E 



^2.26) 



jTo < // I j > 

Third, provided V(x) does not have a tail, there are no cuts, and G(x, y; lj) is real along 
the imaginary axis, which also follows from the reflection symmetry of G. Thus the integral 
along r„i is purely imaginary and does not contribute to G. 
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For the integral along T n 2, the original argument sketched in the previous subsection 
is valid: when Q n is sufficiently large, the integrand goes as exp(— | luj | a), a > 0, and 
therefore vanishes in the limit Q n — > oo; the condition t > t p for case (a) is necessary in 
order that a > 0. 

Next, divide each contour r ra3 into two parts, with | ui | large in one part, and | uoj | small 
in the other. In the first part, the integral also vanishes as the contours expand to infinity, for 
the same reason mentioned previously The complication, as indicated earlier, lies with the 
second part, where | ujj | is not large, and therefore does not provide an exponential damping 
factor. However, in this part the integrand is at most a power of uj and the regulating factor 
e~ WT causes it to vanish as the contours expand to infinity. 

As the contribution from the contour T n vanishes when n — > oo, the Green's function 
G(x,y;t) can be decomposed into a regulated sum of QNM's, namely 

lM)IAii)< ' 7: ' 



G(x, y; t) = Re lim + 



^2.27) 



Thus the prescription is, very simply, to (i) consider only ojr > and (ii) where necessary 
give t a small negative imaginary part r. (The exceptional case where some QNM's lie 
exactly on the imaginary axis is readily handled ||19||.) 

The physical meaning of the regularization is as follows. The Green's function G(x, y; t) is 
a well defined object, independent of the size of the contour V. We attempt to write it as the 
sum of a prompt part arising from the integral on a large semicircle (or other large contour 
linking u = —ioo to uj = +oo), and the sum over QNM's. When the QNM frequencies 
extend to infinity, and there are insufficient powers of 1/u in the integrand, each of the two 
contributions individually diverges as the contour expands to infinity fl7|| . Physically, in 



these circumstances, high frequency QNM's cannot be cleanly separated from the prompt 
contribution, which is hardly surprising. However, the device of letting t — > t — ir effects 
such a clean separation which is practically useful, as illustrated in the numerical examples 
in Sec. IV. 

Mathematically, the regulator e - ^ 7 " can be replaced by any other regulator Ij(r) satis- 
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fying Ij{r) — > 1 as r — > and /j(r) — > as j — > oo. One possible alternative is 

I j (T) = (N-j + l)/N, j<N, 

= 0, j>N, (2.28) 

where iV is the least integer greater than 1/r; this corresponds to defining the divergent 



series by a Cesaro sum [ 20 1 



In practice, one would use a small value of r rather than take r — > + . Consider, for 
example, the regulator e - ^ 7 " and using G to propagate given initial data. Leaving r finite 
incurs an error only for those modes with \ojj\~l/r, and hence only affects the resultant 
wavefunction on length scales below A ~ l/\uj\~r. Thus, so long as only a finite spatial 
resolution is required, not taking r — > + has little effect, as demonstrated in the numerical 
example. 



D. Full line problem 

The proof of completeness can be generalized to include cases with multiple discontinu- 
ities and also to problems on a full line: — oo < x < oo. The latter is physically relevant 
if the background spacetime is due to a Schwarzschild black hole of mass M. Under the 
transformation (1.8), or more explicitly 

x = r + 2Mln(^-l), (2.29) 

the event horizon (r = 2M) maps into x —>■ — oo, while r — > +oo maps to x — > +oo. 
The main idea is still to construct the Green's function G from the two homogeneous so- 
lutions of the KG equation and then examine its high frequency behavior using the WKB 
approximation. Consider a full line problem where there are n discontinuities in V(x) at 
We show only the case a\ < y < x < a n ; the case where the obser- 
vation point is outside (ai,a n ) can be demonstrated similarly. The QNM's are defined by 
(2.1) with the outgoing waves conditions at x — > ±oo. For simplicity we shall assume that 
lim^-too V(x) = 0, and consequently 
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lim fj(x) ~ exp(iu \ x |). 



(2.30) 



In addition, we assume that V(x) has no tail on either side, in the sense of (2.13). 

Introduce two auxiliary functions g_(u),x) and g + {uj,x) which are the solutions 
to the time-independent KG equation Dg± = 0, with the boundary conditions 
linxj^-i-oo [g±(x)e T%UJX ] = 1. As before, the Green's function in frequency space is 



G(x, y; u) = g+(u, x)g-(u, y)/W(u) for y < x, 
= 9+(v,y)g-(u,x)/W(u>) forx<y, 



(2.31) 



where the Wronksian W{uj) = g + (u, x)g'_(u, x) —g' + (ui, x)g_(u, x) is again independent of x. 
Since V(x) has no tail, g±(u, x) are analytic in u>, and hence G is also analytic except at the 
zeros of W(u). If one can show that the integral along the semicircle at infinity vanishes, 
the QNM's will form a complete set, as described in Sec. IIB. 

Let us examine the asymptotic behavior of g±(u,x) along a large semicircle in the half 
plane using the WKB approximation. Assume that a\ < a 2 < ... < a n , then for x > a n 



g+(w,x) ~ exp[il(a n ,x)}g + (u,x = a n ). 



(2.32) 



Now g + (cu, x) in general will consist of two counter-propagating waves for aj < x < aj +1 and 
can be expressed as follows: 



g+(u,x) ~ Ajexp[il(a,j,x)] + Bjexp[— il(aj, x)]. 



The coefficients A~ and B~ can be obtained recursively from the relation 



(2.33) 



1 M 11 (j) M 12 (j)\ 



M 2l {j) M 22 (j) 



x 



(2.34) 



V 



A 



3-1 



exp[i/(aj_i, dj)] 

exp[— i/(aj_i, aj)] 

The boundary condition is A n = g + (cu,x = a n ) and B n = 0. The transfer matrix M(j) 
joining the wavefunction across the discontinuity at x — dj can be obtained by matching 
the boundary conditions at the specified point. At high frequencies 
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M u (j) = M 22 (j) ~ 1, 



(2.35) 



and 



M 12 (j) = M 21 (j)~-R j , 



(2.36) 



where i?j is the reflection coefficient defined by (2.15) and evaluated at x — dj. 

The asymptotic form of g_(u,x) can be similarly expressed in terms of two counter- 
propagating waves for a,j < x < 



g-(u,x) — Cjexp[il(a,j,x)] + Djexp[— il(a,j, x)], 



(2.37) 



where the coefficients C,- and obey the same recursion relation as Aj and Bj, but with 



the initial conditions 



yD 1; 



Mn(l) M 12 , 
M 21 (l) M 22 (l) 



K g-(v,x = ) 



(2.38) 



Then for x < a\ 



g-(u,x) ~ exp[-i/(ai,x)]^_(u;,x = a x ). 



(2.39) 



We now demonstrate the completeness of QNM's for x G a„) as follows. As uj — > oo 
in the lower half plane, the two auxiliary functions can be simplified by keeping only waves 
which grow exponentially in their expansions, and consequently 



g+(u,x) ~ -R n exp[-il(a n ,x)}g + (u,x = a n ), 



(2.40) 



g-(u,x) — —Riexp[+il(ai,x)]g-(u>,x = a\). 



The Green's function is thus given by 



G(x, y; u) ~ exp{i[7(ai, y) + I(x, a n ) - I(a u a n )]}/ (2iu) 



(2.41) 



(2.42) 



for a>i < y < x < a n , and vanishes as uj — > — oo. The completeness of QNM's then follows. 
There is the same problem in the domain \ujj\ = 0(log \wr\), and the same regularization 
needs to be applied. 



20 



III. DIFFERENT SENSES OF COMPLETENESS 



We are now in a position to specify precisely what we mean by completeness, and to 
summarize our results in terms of these definitions. Completeness can be expressed in 
different forms and the meanings may or may not be the same. The first sense (CI) shall 
mean the validity of the following expression, in a certain domain of x, y, t, for the Green's 
function G for problems with an outgoing boundary condition: 

fj(x)f s (y)t '"^iv! 



CI: G(x,y;t) = limRe 



E 



(3.1) 



j>0 - « Ij li » 

In (3.1), the source point is (y,t' = 0) and the observation point is (x,t); for most cases of 
interest in gravitational systems, y < x. The sum is over all QNM's such that the complex 
frequency ujj has Re ujj > 0; the eigenf unctions are fj(x) and <C fj \ fj ^> is the generalized 
norm. The regulating factor Ij could be 



J,(r)=e-^, 



(3.2) 



which corresponds to evaluating G(x,y;t) at a complex t with a small imaginary part. 

We have shown that for discontinuous potentials without a tail, CI holds if (a) < 
y < a < x and t > t p > 0; or (b) < y , x < a and t > 0. We further note that with 
the regulating factor Ij(r), for (x,t) in these two ranges (a) and (b), the sum (3.1) holds 
uniformly and it is valid to differentiate term by term. This leads to 



C2 : G(x, y; t) = Re 



fj^fj^e-^Ijir) 



(3.3) 



b >o < I j J) > 

C2 holds under the same conditions as CI, provided a suitable regulating factor is used. 

The importance of CI and C2 is that they determine the evolution of initial data. For 
t > 0, we have 



<j>(x, t) = Jdy [G(x, y; t)<j>(y, 0) + G(x, y; t)<t>(y, 0) 



(3.4) 



The expansion of G(x, y; t) and G(x, y; t) in terms of QNM's implies that the time develop- 
ment is uniquely determined by the QNM's. Therefore, for discontinuous potentials without 
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tails, if <p(y,0) and <p(y,0) have support only inside (0, a), <p(x < a,t) is completely deter- 
mined for all times by the QNM's, while <p(x > a, t) is completely determined by the QNM's 
for t > t p . 

How is this notion of completeness in terms of evolution related to the notion of com- 
pleteness in terms of the expansion of the delta function? We note that from the defining 
equation and initial condition for G, G(x, y;t = + ) = 5(x — y); so if the domain of validity 
of C2 includes t — > + , then one has in particular 

fM)fj(n)iM) 



C3 : S(x - y) = Re 



,> « h h » 



(3.5) 



which may be seen as the more familiar notion of completeness [[jj 



It is necessary to distinguish these notions of completeness because, unlike the familiar 
case of hermitian systems, completeness in the sense of evolution (CI and C2) is not neces- 
sarily equivalent to completeness in terms of the expansion of the delta function (C3). This 
comes about because each of these holds only in limited regions of spacetime: If y < a < x, 
then CI and C2 only hold for t > t p > 0; as t is bounded away from + , C3 may not be valid. 
Conversely, suppose we know C3 to start with. In the familiar case of hermitian systems, 
such a resolution of the identity operator S(x — y) would allow a decomposition of the initial 
data into a sum of eigenfunctions; attaching phase factors exp(— iujjt) to each of these would 
then give the solution for dynamic evolution. This procedure gives an expression for the 
Green's function G for t > 0, thus proving CI and C2. This is why completeness in the sense 
C3 is usually regarded as the key concept for hermitian systems. However, in the present 
case, because C3 is valid only inside the "cavity" (0, a), the superposition of eigenfunctions 
would only be valid if x G (0, a). Outside this interval, the initial condition would not be 
satisfied by a sum of eigenfunctions and consequently one does not have CI and C2 from 
C3 in general. 

There is another important case for which C3 does not lead to CI and C2, which is 
useful to sketch here. Consider the case V a = ml = linx^oo V(x) > 0. In this case, 



the asymptotic wave number k and the frequency uo are related by k = ycu 2 — The 

22 



asymptotic behavior would be, for example, lim^oo [g(u, x) exp(— ikx)\ = 1. Most of the 
derivation is easily adapted, except that the relation between k and oj leads to an extra 
cut S on the real oj axis. The integral for G(x,y;t) in terms of G is above this cut, so in 
general, there will be this extra cut contribution when the contour is distorted to the lower 
half plane. Thus the QNM sum is not complete in the sense of CI and C2. Nevertheless, 
it turns out that this cut contribution vanishes for t — > + , so the QNM representation of 
S(x — y), i.e., C3 remains valid. 

To see this, define the integral around the cut S 



tional to the product of the regular solutions f(u,x) and f(cu,y). These functions diverge 
exponentially as x or y — > oo along the real axis, since the asymptotic wave number is now 
purely imaginary. However, they are normalizable along the imaginary x or y axis using the 
usual box normalization. It can be shown that G$ is expressible as an integral along the 
imaginary k axis as follows: 




(3.6) 



As both G(x, y\ oj + iO) and G(x, y\ oj — iO) satisfy (2.2), their difference is evidently propor- 




(3.7) 



where f(oj,x) is proportional to f(oj,x) and satisfies the box normalization condition 




(3.8) 



The frequency oj is now given by 



uj = \ m 



2 _ 7,2 1 1/2 



(3.9) 



o 



It is obvious from (3.7) that Gs(x,y;t = + ) vanishes and does not contribute to the QNM 
decomposition of 5(x — y). 
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For open systems as studied here, the importance of C3 is twofold. For initial data 
having support in (0, a), it enables the expansion of initial data in terms of QNM's. It is 
also important in the construction of a time-independent perturbation theory in terms of 
QNM's, an issue we investigated in the case of optical systems |TT[ . 
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IV. NUMERICAL VERIFICATION 



Consider a potential defined on the full line as follows 

V(x) =Vi for < x < a, 

= otherwise, (4.1) 

with V\ > 0. This potential has a step discontinuity (p = 0) and no tail, thus satisfying 
the conditions necessary for our results. Moreover, such a repulsive potential has significant 
leakage, so it is particularly suited for verifying that our results are not merely valid to some 
low order in the amount of leakage. The QNM's of this system are determined by 

fj(x) = exp(iux) for x > a, 

= Aexp(ikx) + Bexp(-ikx) for a > x > 0, 

= -Dexp(— iux) for > x, (4.2) 

where k = \/tu 2 — V±, and the QNM frequency Uj is the j-th solution to 

r(uf exp{2ika) = 1, (4.3) 

with reflection coefficient r(u) = (uj — k)/(u + k) ~ 1/(Auj 2 ) at high frequencies. The result 
(4.3) is a direct generalization of (2.21) by merely replacing R{uj) with r(u) 2 , as waves are 
now reflected successively by the two discontinuities at x = and x — a. The distribution 
of QNM's on the complex uj plane shown in Fig. 3 corresponds exactly to this model, and 
it is easy to show that 

k 3 a ~ J7T - 2t[\og(\j\n/(V 1 1/2 a) + log 2] (4.4) 

for \j\ ^> 1, consistent with (2.22) with R ~ u~ A . Im uj increases with j, implying large 
leakage at high frequencies. 

In the following we shall demonstrate numerically the following approximate equality for 
small r 

25 



Re^ fMmiM ^six-y), (4.5) 
,>o </il/i> 

where iy(r) = exp(— ttyr). We have labeled the QNM's so that j > corresponds to Re 



cjj > 0. Since (4.5) is expected to hold only in a distribution sense [jig) , we integrate 2/ from 
l/i to J/2 







J j( r ) 







Re E - — J% i j J - - yi) - °( x - v*)- ( 4 - 6 ) 

j>0 

Denote the partial sum on the left up to j = J as Sj(x; 2/1, 2/2; r )- Fig. 5a shows \Sj\ versus 
J for a case where x £ (2/1, 2/2), and Fig. 5b shows | Sj — 1 | versus J for x G (2/1,2/2)- In 
each case, curves are shown for several values of r. It is seen that for sufficiently small r (in 
fact r« A = i/2-i/i), the equality indeed holds to a good approximation. These curves 
also show that convergence in j is more rapid for a slightly larger r; this is natural from the 
regulating factor. This property implies that in practice, an extremely small r may not be 
optimal. 

To study the dependence on r, Fig. 6a shows a function y? (x) (e.g., some initial data) 
on the interval (0, a), and Fig. 6b shows the absolute error in representing this function by 
the QNM sum, using r = 10~ 2 , 10~ 3 , 10~ 4 . It is seen that the error converges to zero as 
t — > 0, and that for most purposes, using a small finite r does not matter in practice. The 
QNM sum does not converge if the regulator is removed. 

It is useful to examine through this example why the series needs to be regulated from a 
broader point of view. Let r(u) ~ oj~ q at high frequencies (q = 2 for the present situation). 
By using a similar analysis as that applied in obtaining (2.22), it is straightforward to show 
that Im kjd ~ —q log Next consider the middle region < x < a, and a typical term in 
the product fj(x)fj(y) 

\f j (x)f j (y)\ ~ e ik ^ x+y) ~ ( J - 7r )9(l*/«-i/2|+li//a-V2|+i). (4.7) 

Moreover, it is readily shown that 

« fi I fi » ~ i./^)". (4.8) 
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Since the maximum value of \x/a — 1/2 1 + \y/a — 1/2 1 is 1, the worst behavior of 
\fj(x)fj(y)\/ <C /j | /j 3> is \j\ g . It is seen that the sum (4.5) would not converge, even in 
a distribution sense, without a regulating factor if q > 2. The culprit is Im kjd ~ — glog 
from which one can say unequivocally that the need for regulating the sum is an intrinsic 
property of open KG systems. 

It is worthwhile at this point to reiterate a difference between the wave equation (WE) 
and the KG equation. In each case, if there is a discontinuity in the p-th derivative (of the 
dielectric constant in one case and of the potential in the other case) then 

WE : q = p, 
KG : q = p + 2. 

Thus in the case of the KG equation, the QNM sum requires regulation for all p > 0, whereas 
in the case of the WE, regulation is required only for p > 2. This is again natural when 
seen in the context of (1.10) — a step discontinuity in V is equivalent to a second-order 
discontinuity in the dielectric constant. This difference between the two systems, though of 
a rather technical nature, should nevertheless be kept in mind when using one system to 
draw conclusions about the other. 

The need for a singularity is a feature that might seem, at first sight, to be surprising. 
Our results are valid for a discontinuity however soft (i.e., p however large), but are not valid 
for a C°° potential. Yet there should be little difference between a very soft discontinuity 
and no discontinuity. The resolution of this paradox lies in the regulator — for large p, 
the divergence is more severe, the regulated sum converges much more slowly to the correct 
result as r — > 0, and the use of a small finite r becomes increasingly inaccurate. This 
behavior then connects smoothly to that for a C°° potential, for which the sum may not 
make sense for any r. 
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V. DISCUSSIONS AND CONCLUSION 



In this paper we studied the completeness of QNM's for linearized waves propagating 
in a curved background described by (1.1). It is well known that isolated Schwarzschild 
black holes do not have a complete set of QNM's, while other model systems described by 
(1.1) do have complete sets of quasinormal modes, e.g., the stellar oscillation model of ||. 
The questions we investigate in this paper are: (i) What does completeness mean for the 
QNM's of an open system described by (1.1)? (ii) What properties of the potential would 
characterize those systems with complete sets of QNM's? The answers to these questions 
are spelt out in Sees. II and III. 

We find that there are two important ingredients needed for completeness. First, the 
potential has to have no "tail" in the asymptotic region in the sense of (2.13). Furthermore, 
for potentials behaving asymptotically as 1(1 + l)/x 2 + V(x), where I is an integer, it can 
be shown that the completeness relation holds provided V(x) satisfies (2.13). Details of 
this generalization is given in other context ||21|| . On the other hand, if V(x) decays slower 
than an exponential, the QNM's may not form a complete set. It has been shown that 
for V(x) ~ l/x n with n > 2, the corresponding Green's function G(u) has a branch cut 
along the negative imaginary to axis, and consequently QNM's are incomplete |ID| . The cut 
is related to the late time tail of gravitational waves |ItJ, which decays only as an inverse 
power of time. This in turn verifies that QNM's do not form a complete set. 

Second, the potential has to have discontinuities to provide a demarcation of a finite 
interval, analogous to the boundaries of an optical cavity. The discontinuities can be in any 
finite order of the spatial derivative of the potential. Discontinuities in the potential might 
seem unnatural, however we note that such discontinuities in V commonly exist in many 
models of gravitating systems, e.g., the Price potential as a model for the potential of a 
black hole [p2]| , with the discontinuity representing the peak of the potential; and the stellar 
oscillation model of |§, with the discontinuity representing the stellar surface; as well as 



other stellar models p3|j2^ |. 
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Needless to say there are many models of gravitating systems for which these two condi- 
tions for completeness are not satisfied, with perhaps the most important example being an 
isolated Schwarzschild black hole. For these cases, we note that our study on completeness is 
still relevant in that (i) it provides understanding as to why the quasinormal modes of these 
systems should be incomplete, and (ii) it shows that their QNM's can be made complete by 
a small change in the potential, e.g., by setting the potential V(x) to zero for x > X for a 
large X where V is very small. On the one hand, one does not expect a small change to 
alter the physics much, yet on the other hand, the set of QNM's now becomes complete and 
can be a very useful tool in the analysis of black hole perturbation. This very interesting 
point deserves a more careful discussion, which we now provide. 

Consider a smooth potential, labeled schematically as A; for example this could be the 
Schwarzschild potential. Approximate it by a discontinuous potential B, which vanishes 
rapidly at spatial infinity; this could, for example, be a piecewise constant approximation 
to a Schwarzschild potential P5| . Clearly the latter can be chosen such that the physics 
of the two systems are nearly the same; schematically PHY(A) ~ PHY(B). Yet it is not 
difficult to see, e.g., by reference to the sort of estimates such as (2.22), that the distribution 
of QNM's must be fundamentally different, i.e., QNM(A) ^ QNM(B). Such a qualitative 
difference has been demonstrated numerically in the example of a piecewise constant approx- 
imation to Schwarzschild potential ||25j| . From our point of view, the important difference 
lies not only in the distribution of QNM's, but also in that QNM(A) is not complete, while 
QNM(B) is complete. These differences, both in the asymptotic QNM distribution and in 
the completeness, depend on the order of the discontinuity and not on the magnitude, and 
therefore does not go away even when B is very close to A. This situation is sometimes 
regarded as paradoxical p5| . The paradox is resolved if one realizes that the QNM's are not 
complete for system A, and their sum does not correctly describe the physics, i.e., PHY(A) 
^ QNM(A). Of course, the results in this paper implies PHY(B) = QNM(B). Thus, in all, 
we have QNM(A) ^ PHY(A) ~ PHY(B) = QNM(B). 

This remark opens up a particularly intriguing possibility. If we want to approximate 
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the time development of A by a discrete sum of QNM's, we should not use the QNM's of 
the system itself. Instead, the QNM's of B, which are complete, can describe the evolution 
of A to a good approximation. Investigations along this line will be reported elsewhere. 
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FIGURE CAPTIONS 



Fig. 1. The integrand in (2.10) is defined by the contour L = Li + L 2 , where L x = (0, a) 
and L 2 is the contour shown in this figure. 

Fig. 2. Contributions from the poles, the semicircle, and the singularities of / and g (which 
in general form a cut on the negative imaginary axis). 

Fig. 3. Distribution of QNM's for the model in Sec. IV with a = 1 and V\ = 100. 

Fig. 4. Contributions from the poles, and the rectangular contour r„ = T„i U T n2 U r n3 . 

Fig. 5. (a) Plot of log \Sj\ vs J for a case where x £ (yi,^)- (b) Plot of log \Sj — 1| vs J 
for a case where x G (yi,y2)- In each case a — 1, V± — 100, J is in thousands, and Sj 
is the partial sum in (4.6), evaluated for r = 5 x 10 -3 , 5 x 10~ 4 , 2 x 10 -4 for lines 1, 
2, 3 respectively. The partial sums are in fact fluctuating functions of J, and the lines 
shown are smooth envelopes representing upper bounds. These curves show that (i) 
at fixed r, the sum over j converges; (ii) the resultant error vanishes as r — > 0. 

Fig. 6. (a) A function (p Q (x) = sin(7re/a) defined on (0,a). (b) The absolute error in 
representing this function by the QNM sum, using r = 10~ 2 (circles), 10~ 3 (squares), 
10~ 4 (triangles). In this example, a = 1 and V\ = 100. 
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